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Some aspects, both experimental and theoretical, of extracting the neutron skin AR from 
properties of isovector giant resonances are discussed. Existing proposals are critically reviewed. 
The method relying on the energy difference between the GTR and IAS is shown to lack sensitivity 
to AR. A simple explanation of the linear relation between the symmetry energy and the neutron 
skin is also presented. 

1 Introduction 



>■ | Accurate experimental data on the neutron skin in neutron rich nuclei would allow to further constrain 
model parameters involved in the calculations of the nuclear symmetry energy [lj. The latter plays 
a central role in a variety of nuclear phenomena. The value rs 30 MeV of the nuclear symmetry 
energy S(p ) = a 4 + ^(p — p ) + ... at nuclear saturation density p ~ 0.17 fm~ 3 seems reasonably well 
established. On the other hand, the density dependence of the symmetry energy can vary substantially 



t — ! with the many-body approximations employed. 
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Several authors have pointed out [21 [3] a strong correlation between the neutron skin, AR = yj (r 2 ) n — 



(r 2 ) p = R n —Rp, and the symmetry energy of neutron matter near saturation density. In the framework 
of a mean field approach Furnstahl [3] demonstrated that in heavy nuclei there exists an almost linear 
empirical correlation between theoretical predictions in terms of various mean field approaches to S(p) 
(i.e., a bulk property) and the neutron skin, AR (a property of finite nuclei). 

This observation has contributed to a renewed interest in an accurate determination of the neutron 
skin in neutron rich nuclei. Besides, a precise value of the neutron skin is required as an input in several 
processes of physical interest, e.g. the analysis of energy shifts in deeply bound pionic atoms 0], and 
in the analysis of atomic parity violation experiments (weak charge) [5]. It is worth to stress that to 
experimentally determine the skin in heavy nuclei is extremely challenging as AR is just about few 
percents of the nuclear radius. 

The present contribution is partially based upon the results published previously in [BJ. 



2 Relationship between the symmetry energy and AR 

Brown [2] and Furnstahl [3] have pointed out that within the framework of mean field models there 
exists an almost linear empirical correlation between theoretical predictions for both 04 and its density 
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dependence, po, and the neutron skin AR in heavy nuclei. This observation suggests an intriguing 
relationship between a bulk property of infinite nuclear matter and a surface property of finite systems. 
Here, following the analysis of [6J, this question is addressed from a point of view of the Landau-Migdal 
approach. 

Let us consider a simple mean-field model with the Hamiltonian consisting of the single-particle 
mean field part Hq and the residual particle-hole interaction H p -h'- 

H = H + H P - h , H ph =E(F' + G'a a a b )r a f b 5(f a - r b ), (1) 

a>b 

H = E(T a + U(x a )), U(x) = U (x) + U^x) + U c (x), (2) 

a 

U (x) = U (r) + U 9O (x); U^x) = \S vot {r)r^- U c {x) = \U c {r){l - r< 3 )). (3) 

Here, Uq(x) is the phenomenological isoscalar part of the mean field potential U(x) (x = {r,a, r}), 
Uo(r) and U so (x) are the central and spin-orbit parts, respectively; F' and G' are the phenomenological 
Landau-Migdal parameters. The isovector part U\(x) and the Coulomb mean field Uc(x) are both 
calculated consistently in the Hartree approximation, S pot (r) is the symmetry potential (r-dependent 
symmetry energy in finite nuclei). 

The model Hamiltonian H (CQ) preserves the isospin symmetry within the RPA if a selfconsistency 
relation between the symmetry potential and the Landau-Migdal parameter F' is fulfilled: 

S pot (r) = 2F'n(-)(r), (4) 

where n^(r) = n n {r) — n p (r) is the neutron excess density. Thus, in this model the depth of the sym- 
metry potential is controlled by the Landau-Migdal parameter F 1 (analogous role plays the parameter 
g 2 p in relativistic mean field models). S^otM is obtained from Eq.(|4]) by an iterative procedure; the 
resulting dependence of AR on the dimensionless parameter /' = F'/(300 MeVfm 3 ) shown in fig. CD 
indeed illustrates that AR depends almost linearly on /'. Then with the use of the Migdal relation 
04 = ^(1 + 2f) [7] relating the symmetry energy and /', a similar, almost linear, correlation between 
04 and AR is obtained. 

To get more insight in the role of /' we consider small variations 5F'. Neglecting the varia- 
tion of n^~\r) with respect to 5F', the corresponding linear variation of the symmetry potential 
is 5S pot (r) = 25F'nS \r). Then in the first order perturbation theory, such a variation of S pot 
causes the following variation of the ground-state wave function \50) = 5F' ^Eo-e^ with V 

labeling the eigenstates of the nuclear Hamiltonian and a single-particle operator defined as 

N(~i = n^\r a )r^ 3 \ Consequently, the variation of the expectation value (0|V^~)|0) = NR^ — ZRl 

a 

of another single-particle operator V^' = r a r a 3 ' ) can be written as 

a 

R P 5{AR) = 5F ■ Eq-E s ' ^ 

In practice the sum in Eq. (JHJ) is exhausted mainly by the isovector monopole resonance (IMR) which 
high excitation energy (about 24 MeV in 208 Pb) justifies the perturbative consideration. Eq. (jHJ) is able 
to reproduce directly calculated S(AR) shown in Fig. [T]with the accuracy of about 10%. As a result, 
a simple microscopic interpretation of the linear correlation between AR and Landau parameter F' is 
obtained. 
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Figure 1: Neutron skin in 208 Pb versus the Landau- Migdal parameter /'. 



3 Extracting neutron skin from properties of isovector giant 
resonances 

Parity violating electron scattering off nuclei is probably the least model dependent approach to probe 
the neutron distribution [8]. The weak electron-nucleus potential is V(r) = V(r) + 75,4(7-), where 
the axial potential A(r) = i02Pw{ r )- The weak charge is mainly determined by neutrons pw( r ) — 
(1 — 4 sin 2 9w)Pp(r) — p n (r), with sin 9w « 0.23. In a scattering experiment using polarized electrons 
one can determine the cross section asymmetry [5] which comes from the interference between the A and 
V contributions. Using the measured neutron form factor at small finite value of Q 2 and the existing 
information on the charge distribution one can uniquely extract the neutron skin. Some slight model 
dependence comes from the need to assume a certain radial dependence for the neutron density, to 
extract R n from a finite Q 2 form factor. 

However, the best claimed accuracy of the experimental determination of neutron radii would be on 
the level of 1%, that translates to relatively large uncertainty of 20-30% in the neutron skin. On such 
accuracy level, some indirect experimental probes of AR still can be competitive. 

A variety of experimental approaches have been employed to obtain indirect information on AR. 
To some extent all the analysis contain a certain model dependence, which in many cases is difficult 
to estimate quantitatively. For choosing an indirect probe it is very important to address the question 
how sensitive is the proposed physical quantity with respect to a variation of AR in a single nucleus. 
The higher is the sensitivity, the better is the choice of the correlation for the indirect deducing AR 
from the measured values. 

It is not intended here to give a comprehensive review of the existing methods. In particular, the 
results from the analysis of the antiprotonic atoms, elastic proton and neutron scattering reactions, and 
the pygmy dipole resonance are completely left out. Here, special emphasis will be put on proposals to 
provide accurate information on the neutron skin from properties of isovector giant resonances. 
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3.1 Spin-dipole Giant Resonance 



In [9] it has been proposed to utilize the excitation probability of the spin-dipole resonance in charge 
exchange reactions for determining the neutron skin. The method has been applied to obtain information 
on the variation of the neutron skin in the Sn isotopes [9j . For the relevant operator, J2a T t [&a ® 
fJjM, (J = 0, 1, 2) the summed AL = 1 strength is 

S(~) - £(+) = C(NR 2 n - ZRl). (6) 

Here and are the spin-dipole total strengths in (3^ and (3^ channels, respectively; C is the 
factor depending on the normalization of the spin-dipole operator (in the definition of Ref. [2] C = 1/47T, 
we use here C — 1). Because could not be measured experimentally, the model-dependent energy- 
weighted sum rule was invoked in the analysis of [9 J to eliminate S^ + \ However, the used analytical 
representation for the sum rule was oversimplified and led in some cases, e.g. for 208 Pb, to absurdly 
negative S^ + \ In [TU] another way was proposed, namely, to use for the analysis the ratio S^/ 
calculated within the pn-RPA. The parameterization of the RPA calculation results for tin isotopes in 
the form 

S-M/S-R = o.388 - 0.012(iV - Z) 

was used later in [TT] to reanalyze the experimental data and led to a marked change in the extracted 
ARs. 

Let us now assess the experimental accuracy for needed to determine the neutron skin to a 
given accuracy. Putting = (that seems to be a very good approximation for 208 Pb) and one has 

S-H = (N - Z)R 2 p + 2NR P AR. (7) 

The ratio of the second term on the rhs to the first one in case of 208 Pb is 

2NAR/((N - Z)Rp) w 5.7AR/R P . 

Therefore, for R p = 5.5 fm and AR = 0.2 fm the second term is only 25% of the first one and one needs 
5% accuracy in to determine AR with 20% accuracy. Because the SD strength is spread out and 
probably has a considerable strength at low-energy, the results for the AR can be only considered as 
qualitative with a relatively large uncertainty (up to 30-50%). 



3.2 Isobaric analogue state 

The dominant contribution to the energy weighted sum rule (EWSR) for Fermi excitations by the 
operator = J2 a T a comes from the Coulomb mean field 

(EWSR) F = J U c (r)n(-\r)d 3 r, (8) 

The Coulomb mean field Uc{r) resembles very much that of the uniformly charged sphere, being inside 

Ze 2 

a nucleus a quadratic function: Uc(r) = 77rr(3 — (r/R c ) 2 ), r < R c . It turns out that if one extends 

2R C 

such a quadratic dependence also to the outer region r > R c (instead of proportionality to R c /r), it 
gives numerically just a very small deviation in (EWSR)f (less than 0.5%, due to the fact, that both 
the difference and its first derivative go to zero at r = R c and n'~'(r) is exponentially decreasing for 
r > R c ). Using such an approximation, one gets: 

WSV^IN-Q^l-^^L-) (9) 
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with Ac = - , and ' given in Eq.([7|). 

2R C 

Since the IAS exhausts almost 100% of the NEWSR and EWSR, one may hope to extract from 
the IAS energy. However, the term depending on contributes only about 20% to (EWSR)f, and as 
a result, the part of depending on AR contributes only about 4% to (EWSR)f (in 208 Pb). ^From 
the experimental side, the IAS energy can be determined with unprecendently high accuracy, better 
than 0.1%. Also, from the experimentally known charge density distribution the Coulomb mean field 
Uc{r) can be calculated rather accurately, and hence one can determine the small difference between 
Eqs.flHD an d @- But at the level of 1% accuracy several theoretical effects discarded in Eq.([H]) come 
into play that makes such an accurate description of the IAS energy very difficult (the Nolen-Schiffer 
anomaly) . 

Also in [12J it was stated that the Coulomb displacement energies (CDE) are sensitive to AR. A 
gross estimate AR = 0.80(5) (N — Z)/A fm was obtained from a four-parameter fit of the experimental 
R p and observed mirror CDE's. The authors claimed 127 keV to be the rms error of the fit, but they 
assumed the nuclear wave functions calculated within the Nuclear Shell Model to be isospin pure. Thus, 
the important effect of the Coulomb mixing of the IAS and the IMR was not taken into account, which 
is known to decrease the IAS energy by a few percents. Therefore, the Nolen-Schiffer anomaly does not 
seem to have been resolved yet. 

3.3 Is the energy spacing between GTR and IAS a good candidate for 
determining the neutron skin in isotopic chains? 

In a recent paper [13] a proposal has been put forward to use the isotopic dependence of the energy 
spacing, AE, between the Gamow- Teller resonance (GTR) and the IAS as a tool for determining the 
evolution of the neutron skin in nuclei along an isotopic chain. Here, we would like to present some 
physical arguments which question the physical relevance of this method. 

The authors of [13] have used the fact that both functions, AR and AE, are monotonic functions 
(increasing and decreasing, respectively) of the neutron excess (N—Z) to state that "isotopic dependence 
of the energy spacings between the GTR and IAS provides direct information on the evolution of neutron 
skin-thickness along the Sn isotopic chain" . Arguing in such a way one can find a correlation between 
any two monotonic functions of a single physical parameter and plot them as a function of one another 
like is done in Fig. 2 of [13] 0. However, it does not imply automatically a real physical correlation 
between the functions which are determined also by many other model parameters which are kept fixed 
while performing calculations (the calculations in [TJ] have been performed within the relativistic mean 
field (RMF) and relativistic QRPA (RQRPA) approaches). 

Again, the relevant question to be addressed is how sensitive is one physical quantity with respect to 
a variation of another in a single nucleus ? In other words, one has to evaluate what variation of AE 
is produced by varying AR in a single nucleus. Imaging an extreme situation (which is actually not 
far from reality) that AE were not sensitive to AR at all, one would get by varying AR a family of 
different calculated dependences (like shown in the upper panel of Fig. 2 of [13]) which would give no 
clue about the real dependence seen in nature. 

Thus, it is quite important to understand the physical reasons which cause the energy splitting 
between the GTR and the IAS. It is well-known that if the nuclear Hamiltonian possessed Wigner 
SU(4) symmetry then the GTR and the IAS would be degenerate, AE = 0. In such a case any 
variation of AR, not violating the symmetry, would not affect AE at all. However, it is also known 
that the spin-isospin SU(4) symmetry is broken in nuclei. Hence, AE is determined by those terms in 

^ote that, to avoid confusion in comparing the measured and calculated dependences, the authors should have plotted 
the experimental points in the upper panel as the function of the measured AR rather than calculated AR and should 
have added the horizontal error bars to them reflecting the experimental uncertainty in AR shown in the lower panel. 
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the nuclear Hamiltonian which violate the symmetry. Their qualitative and semi- quantitative estimates 
in terms of the energy weighted sum rules for the Gamow- Teller (EWSRqt) an d the Fermi (EWSRp) 
excitations have been already known for more than 20 years (see, e.g., [33]). The analysis of these 
authors as well as a quantitative analysis performed recently in [15] has shown that there are three 
basic sources in the Hamiltonian which violate SU(4) symmetry and contribute to the difference of the 
sum rules: spin-orbit mean field and both particle-particle and particle-hole residual charge-exchange 
interactions. One sees that none of the sources explicitly refers to the symmetry potential, to which 
AR is especially sensitive. 

EWSR G t - EWSR F 

An estimate of AE as AE = can be calculated according to ITS! in the 

Sn isotopes. From the sources violating SU(4) symmetry, spin-orbit mean field represents the major 
one and contributes about 5 MeV to the splitting. The contribution of the particle-hole interaction is 
negative and about 1-2 MeV in the absolute value. The contribution of the particle-particle interaction 
is rather difficult to evaluate (due to uncertainty in the strength of the spin-dependent particle-particle 
interaction) but it seems to be of minor importance (very probably no more than 0.5 MeV, especially 
for large (N — Z)) and can safely be neglected. 

Now let us turn to the discussion of the sensitivity of the contributions to the variation of AR. We 
could reproduce the corresponding analytical expressions from [15] explicitly, but it is enough for our 
purpose just to mention that the dominating contribution to AE from the spin-orbit mean field is given 
by its expectation value in the ground state and is determined basically only by the unfilled spin-orbit 
doublets. This expectation value is completely insensitive to the variation of AR. 

Within the Landau-Migdal approach described above, the particle-hole contribution 

AE ph = 2{G '~P l{^\r)fSr (10) 



N-Z 

is given by the product of the volume integral of the neutron excess density squared and the difference 
of the p-h strengths G' and F' [15]. In the S , f/(4)-symmetric limit one has G' = F' and AE p h = 
explicitly. Still, in this limit one has a freedom to choose different F' that produces a variation in AR, 
similar to shown in fig 1. Therefore, as already mentioned, one can get no clue about the actual AR 
from AE = in the £77(4)-symmetric limit. 

In a realistic situation G' ^ F' (/' = 1.0 and g' = 0.8 were taken in [UJ), but AE p h depends only on 
the difference G' — F'. One usually fixes G' in order to reproduce the GTR energy in some nuclei (the 
authors of [13] have followed this way, too) and possible information from AE about the absolute value 
of F' is lost. Furthermore, one can a priori think that a degree of violation of the SU(4) symmetry 
should be a sort of a fundamental property of the residual interaction. Therefore, the difference G' — F' 
should stay more stable in different models as compared to some possible variation of F' producing 
different AR. 

Considering value of G' — F' fixed, one can employ a simple model varying only p n {r) to see how a 
change of AR affects AE p h via variation of the neutron excess density n^~\r). A small variation of p n (r) 

can be approximately represented as 5p n (r) = — rr"(3p n ( r ) + R ^7 - )> where 8R n is a change of the 

R n dr 

rms neutron radius R n , R is the nuclear radius (with R 2 n « 0.6i? 2 ). Assuming the proton and neutron 

, •• , • • , 1 1 r- , • • SAE ph JRn„ N + Z-2-/N, , 

densities be constant inside a nucleus, the final estimate is — — — = —3 (H , where 

AE ph R n y N-Z h 

7 = 7iH(i2)/ n (-)(0). 

Thus, in Sn isotopes with the experimental charge radii about R p =4.6 fm a rather significant 
variation of AR about 0.1 fm, that is of the order of magnitude of 5R, would cause A£ / fe = 0.3 and 

0.15 for 112 Sn and 132 Sn, respectively (7 = 0.5), that corresponds to the absolute change about 



AE 



P h 



0.3 MeV in AE, to be compared with the experimental uncertainties in AE of the same order. It 
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is clear that to draw any conclusion about AR from the measured AE would be premature. Even if 
the experimental errors in AE were exactly zero, the accuracy of the theoretical model itself would be 
hardly believed to be of the necessary level. For instance, apart from the obvious uncertainties in the 
isotopic dependence of the spin-orbit potential, the GTR does not exhaust 100% of the corresponding 
sum rules and the shell-structure effects such as configurational and isospin splitting of the GTR can 
have some effect on the calculated GTR energy. 

It is also noteworthy that, in spite of the claimed self-consistency of the calculations, the slope of the 
calculated isotopic dependence of the IAS energy is about 3 times larger than the experimental one (see 
inset in Fig. 1 of p3]). Note, that the isospin self-consistent continuum- QRPA calculations of [15] were 
able to nicely reproduce the slope (while overall underestimated the IAS energy by about 0.5 MeV, the 
well-known Nolen-Schiffer anomaly). 

To conclude, we believe that the suggested in [13] method to deduce the neutron skin from the 
energy spacing between GTR and IAS is rather questionable in its origin and does not fairly provide 
"direct information on the evolution of neutron skin-thickness" . 



4 Some implications of AR 

In several processes of physical interest knowledge of AR plays a crucial role and in fact a more accurate 
value could lead to more stringent tests: 

(i) The pion polarization operator [4] (the s-wave optical potential) in a heavy nucleus II(u;, p p , p n ) = 
—T + (u))p—T~(uj)(p n —p p ) has mainly an isovector character (T + (m 7r ) ~ 0). Parameterizing the densities 
by Fermi shapes for the case of 208 Pb the main nuclear model dependence in the analysis comes from 
the uncertainty in the value of AR multiplying T~ . 

(ii) The parity violation in atoms is dominated by Z— boson exchange between the electrons and the 
neutrons [5j. Taking the proton distribution as a reference there is a small so-called neutron skin (ns) 
correction to the parity non- conserving amplitude, 5Ep^ c , for, say, a 6sy 2 — > 7s 1/2 transition, which is 
related to AR as (independent of the electronic structure) 



Enwc 7 R 



--{aZY — . (11) 



>pnc 

In 133 Cs it amounts to a 5E/E w —(0.1 — 0.4)% depending on whether the non-relativistic or relativistic 
estimates for AR are used [5]. The corresponding uncertainty in the weak charge Qw is —(0.2 — 0.8)c. 

(iii) The pressure in a neutron star matter can be expressed as in terms of symmetry energy and its 
density dependence 

P(p, x) = p^^^l = p\E\p, 1/2) + S'(p)(l - 2xf + ...]. (12) 

By using beta equilibrium in a neutron star, p e = p n — p p = — , and the result for the electron 
chemical potential, p e = 3 / 'Ahcx{2>-K 2 px) 1 ^ 3 , one finds the proton fraction at saturation density, po, to 
be quite small, xq ~ 0.04. Hence, the pressure at saturation density can be approximated as 

P(p ) = p s (l - 2x )(p S'(p )(l - 2x ) + S( Po )x ) ~ plS'(po). (13) 

At higher densities the proton fraction increases; this increase is more rapid in case of larger p [1] . While 
for the pressure at higher densities contributions from other nuclear quantities like compressibility will 
play a role in it was argued that that there is a correlation of the neutron star radius and the pressure 
which does not depend on the EoS at the highest densities. Numerically the correlation can be expressed 
in the form of a power law, R M ~ C(p, M)( Me ^_ a ) a25 km, where C(p = 1.5p 8 ,M = lAM solar ) ~ 7. 
This shows that a determination of a neutron star radius would provide some constraint on the symmetry 
properties of nuclear matter. 
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5 Conclusion 



In this contribution we discuss some aspects of extracting the neutron skin from properties of isovector 
giant resonances and critically review existing proposals. The theoretical method relying on the energy 
difference between the GTR and IAS is shown to lack sensitivity to AR. It is also shown that the phe- 
nomenological, almost linear, relationship between the symmetry energy and the neutron skin in finite 
nuclei, observed in mean field calculations, can be understood in terms the Landau-Migdal approach. 
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